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To seek a deeper understanding of wave propagation on tlie Schwarzschild spacetime, we investi- 
gate the relationship between (i) the lightcone of an event and its caustics (self-intersections), (ii) 
the large-/ asymptotics of quasinormal (QN) modes, and (iii) the singular structure of the retarded 
Green function (GF) for the scalar field. First, we recall that the GF has a (partial) representation 
as a sum over QN modes. Next, we extend a recently-developed expansion method to obtain asymp- 
totic expressions for QN wavefunctions and their residues. We employ these asymptotics to show 
(approximately) that the QN mode sum is singular on the lightcone, and to obtain approximations 
for the GF which are valid close to the lightcone. These approximations confirm a little-known pre- 
diction: the singular part of the GF undergoes a transition each time the lightcone passes through 
a caustic, following a repeating four-fold sequence. We conclude with a discussion of implications 
and extensions of this work. 



I. INTRODUCTION 

Characteristic damped resonances of black holes 
known as quasinormal ( QN) modes have received much 
attention over the last three decades pillSj. QN modes 
are expected to play a key role in strong-field processes 
in black hole dynamics (see e.g. [TBl - fI5] ). For example, 
after the merger of two black holes, the composite system 
undergoes a 'ringdown' phase during which it sheds its 
asymmetries through gravitational radiation. This phase 
is well-described in the linearized approximation as a 
black hole radiating via the least-damped, low-multipole 
QN modes. 

A black hole QN mode corresponds to a complex fre- 
quency, ujQ = 27r(l/r — i/r), where T is the period of 
oscillation and r is the decay timescale. The frequencies 
of the QN spectrum are labelled by integers, e.g. the mul- 
tipole I and overtone n > (in Schwarzschild spacetime). 
The QN spectrum depends only on the properties of the 
field (e.g. its spin s) and the black hole geometry (e.g. 
mass M, angular momentum J), i.e. w;„(s, J, M). Hence, 
the (future) observation of QN modes in gravitational- 
wave data would allow for the direct inference of mass 
and angular momentum of black holes [17' . 

Most studies of QN modes have focussed on the chal- 
lenge of determining the QN frequency spectrum. Rather 
fewer studies have tackled the more demanding question 
of how (and how much) QN modes are excited by a given 
initial data (see notable exceptions [20H25] ). Although 
the QN formalism is well understood [20], such calcula- 
tions face certain technical challenges. For example, un- 
like normal modes in conservative systems, QN modes do 
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not form a complete set (i.e. they provide only a partial 
description of the system). Furthermore, QN mode sums 
may not be strictly convergent. The aim of this work is to 
show that, despite technical issues, the QN mode formal- 
ism can be used to gain geometric insight into intriguing 
features of wave propagation on curved spacetimes. In 
this study, we consider the simplest case: a non-rotating 
(J = 0) black hole of mass M in the Schwarzschild coor- 
dinate system {t, r, 9, (f)} described by line element 

ds^ = -f{r)dt^ + f-\r)dr^ + r^ {dO^ + sm^ ddc/)^) , (1) 

where f{r) = 1 — 2M/r. The event horizon is at r = 2M. 

To motivate this study, let us begin with a simple 
thought experiment, in which a 'beacon' at spatial coor- 
dinate ~a^', close to a black hole, emits an instantaneous 
flash of light at t = 0. A moment later, a,tt = 6t, the light 
forms an almost-spherical wavefront around it' . As t in- 
creases, the wavefront propagates outwards, describing a 
'lightcone': a 3D null hypersurface embedded within the 
4D spacetime. Locally, each part of the wavefront prop- 
agates along a null geodesic of the spacetime. The pres- 
ence of the black hole warps the wavefront, and the part 
of the wavefront that passes close to the 'unstable photon 
orbit' at r = 3M may orbit the black hole multiple times. 
This observation has several concomitants: (i) the light- 
cone intersects itself along null lines in spacetime called 
caustics; the flrst caustic forms after a time At, where 
At is the shortest time it takes for a light-ray emanating 
from it' to pass to the diametrically opposed side of the 
black hole; (ii) the wavefront passes through the spatial 
point at it' (or indeed any point) an infinite number of 
times, at t = ^1,^21 ■ • • (where ti < t2 < ■ ■ .), although 
each 'echo' will be fainter than the last; and (iii) a black 
hole acts as a very strong gravitational lens, generating 
an infinite number of images of the same source [26H28] . 
Some aspects of this thought-experiment are illustrated 
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in Fig.[T} 

In this paper, we study the retarded Green function 
which provides a mathematical description of wave prop- 
agation. In flat spacetime, the retarded Green function 
(GF) has support only on the lightcone itself. In curved 
spacetimes, it is well-known that the GF also has (non- 
singular) support within the lightcone, due to backscat- 
tering from the curved geometry [301 HI] ■ It is perhaps 
less well-known that the singular part of the GF is mod- 
ified when caustics are encountered The princi- 
pal aim of this paper is to show that the singular part 
of the GF on the Schwarzschild spacetime may be well- 
understood through the large-Z asymptotics of the QN 
mode representation. 

In Ref. 'TT (henceforth Paper I) we developed an ex- 
pansion method to investigate the large-/ asymptotics 
of the QN frequency spectrum of spherically-symmetric 
black holes. We showed that the frequencies of the 
lightly-damped modes of the Schwarzschild black hole are 
given by the expansion 
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where /3 = 1 — and 



L = 1 + 1/2, TV = 71 + 1/2, 
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and / > |s| and n > are the angular momentum and 
overtone numbers, respectively. Higher-order terms in 
the Schwarzschild expansion were given in Paper I. Ex- 
pansions for some other spacetimes were computed in 

[35M]. 

In this work, we develop the expansion method fur- 
ther to investigate the asymptotics of the so-called 'exci- 
tation factors' Bin (defined in Sec. Ill I, which provide a 
stepping-stone to the properties of the Green function. 

The remainder of this paper is organised as follows. 
In Sec. [n] we outline the basics of black hole perturba- 
tion theory and quasinormal modes. In Sec. Ill we recap 



the theory that leads to a (partial) decomposition of the 
Green function in a series of QN modes. In Sec. |IV| we 
derive some key results for the large-/ asymptotics of QN 
wavefunctions and excitation factors. In Sec. |V| we em- 
ploy the asymptotics to investigate the singular structure 
of the Green function in the Schwarzschild spacetime. We 



conclude in Sec. VI with a discussion of the applications 
of the method. 



II. BASICS 

The evolution of scalar-field, electromagnetic and grav- 
itational perturbations on the Schwarzschild spacetime 
are governed by wave equations [321 EHI ■ The simplest 
example is the (minimally-coupled) scalar field in the 



vacuum exterior of a black hole, which satisfies a curved 
space Klein-Gordon equation: 
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da 
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Here g = —r'^ src? 9 is the metric determinant, and g'^" 
is the contravariant version of the metric. The symme- 
try of the spacetime allows for a complete separation of 
variables, i.e. 



$(i, = ^'Y^ Yim{9, (fycim J ui^{r)e'' 



*doj, (5) 



where Yim{d, (t>) are spherical harmonics. The radial func- 
tion uicj (r) satisfies the second-order ordinary differential 
equation 

" 2 ^ fKl + '^) 2A//3V 



.W=0, (6) 



where, as above, /(r) = 1 — 2Af/r, and, in addition, we 
have introduced the standard 'tortoise coordinate' = 
r-l-2Mln(r/2M-l). Here /3 = l-s^ with \s\ = 0,1, 2 for 
scalar, electromagnetic and gravitational perturbations 
[5], respectively. On a black-hole spacetime, physically 
valid solutions are purely 'ingoing' at the event horizon. 
Let us define 'ingoing' solutions u\^{r) via 
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where AiJ and A^^^ are complex constants. Quasinor- 
mal modes (QNM) are modes with complex frequencies 
LUin such that A^J^ = (where n = 0, 1, ... is the over- 
tone number). In other words, QNMs are purely ingoing 
at the horizon and purely outgoing at spatial infinity. 

The spectrum of QNMs of the Schwarzschild space- 
time is illustrated in Fig. [2j The plot highlights some 
key properties. Firstly, all QNM frequencies have a neg- 
ative imaginary part (in other words, the amplitude of 
all QNMs decay with time). Secondly, the frequencies 
are approximately evenly-spaced in / and n. In the limit 
n — > oo (for fixed I) the spacing uji^n+i ~ '^i,n approaches 
an imaginary (and /-independent) constant. In the limit 
/ — > oo (for fixed n) it is well-established [3] that the spec- 
trum has the asymptotic form aJi„ w LSI — iTVA, where 
L and N are defined in Eq. ([s]), while VL = {^/27M)-^ is 
the orbital frequency on the 'photon sphere' at r = 3Af , 
and A = {^/27M)~^ is the so-called Lyapunov exponent 
(i.e. decay timescale) [38| for orbits in the vicinity of 
this sphere. In essence, the asymptotic properties of the 
QNM spectrum are linked to the properties of the pho- 
ton sphere; this observation extends also to more general 
black hole spacetimes [15^. In Paper I, we developed an 
asymptotic expansion method to explore this link math- 
ematically, and applied the method to obtain high-order 
expansions in inverse powers of L for the low-overtone 
modes of the spectrum, e.g. Eq. In this work, we 
will extend the method to investigate wave propagation 
in the large-L limit. 



x/M 



FIG. 1. Left: Light rays orbiting a Schwarzschild black hole. This plot shows that two points at fixed spatial coordinates in 
the vicinity of a black hole are linked by multiple null geodesies (i.e. light rays). Here we show the 'direct' ray (ray 0) [passing 
through = it/2] and two orbiting rays [passing through A(p = ?)Ti/2 (ray 1) and A<j!> = 57r/2 (ray 2)] in the xy plane. 
Right: The self-intersecting light cone (reproduced from Fig. 1.1 of |29| : see also [26]). The right plot shows the lightcone of a 
spacetime point (at the apex of the cone) in the vicinity of a black hole (the horizon of the BH is visible as the smaller circle) . 
The coordinate time t runs vertically downwards, and one spatial dimension (z) has been suppressed. Note the formation of 
caustics (lines along which the lightcone intersects itself). See also Fig. [s] and [t] 
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FIG. 2. The gravitational quasinormal mode spectrum of the 
Schwarzschild black hole. The plot shows the real and imagi- 
nary parts of the gravitational (|s| =2) QNM frequencies, for 
angular momenta i = 2 ... 10 and overtones n = ... 11. 



III. THE GREEN FUNCTION AND 
QUASINORMAL MODE SUMS 

The role of QNMs in wave propagation may be ap- 
preciated by examining the contribution of QNMs to the 
retarded Green function [2Ql [22] . Here, for simpHcity, we 
will consider the retarded Green function for the scalar 
field on Schwarzschild spacetime, Grot(a^, x'), which is de- 



OxG,:otix, x) = (5'*(X - x) 



(8) 



(where S"^ is a covariant version of the Dirac delta 
distribution) with appropriate causal conditions. The 
Green function may be used to solve two types of prob- 
lem. Firstly, if a field <^{x) is generated by a source 
S{x) of compact support, i.e. □$(a;) = S{x), then 
the field is found from the spacetime integral <I>(a:) = 
/ G^ct{x^ x')S{x')d'^x' . Secondly, if an unsourced field sat- 
isfies a Cauchy initial value problem then the field may 
be found via an integral over a hypersurface. For exam- 
ple, if at i = the field (in the exterior region of the BH) 
is given by $|(=o = ^oi^) and dt^\t=o ~ then the 
field at later times t is given by 



(9) 



Ht,x)^ j [G,.et(t,a3;0,a;')3'oK) + 
dtG,,t{t.x-{),x')'^Q{x')](i'x' 



The retarded Green function on Schwarzschild space- 
time may be expressed in terms of an inverse Laplace 
transform. 



^ oo 

{x,x')=- ^^(2; + l)P,(cos7); 



+oo-|-ic 



Gz^(r,r')e-^"(*-*')da;. (10) 



oo+ic 



Here x,x' are spacetime points at radii r, r', separated 
by coordinate time t — t' and spatial angle 7, [where 



4 



(2) 



o o o o 



(3) 



'(1) 



FIG. 3. Integration contour in the complex frequency plane. 
In Eq. (10 1, the GF is expressed as an integral in the upper 



half of the complex-frequency plane. In principle, for T > 
[Eq. [15] , this may be computed by closing the contour in the 
lower half-plane. This gives three contributions: (1) a high- 
frequency arc, (2) a sum over the residues of the QN modes 
(which are poles of the GF, indicated by the dots above), and 
(3) an integral along the branch cut which runs along the 
negative imaginary-frequency axis. 



cos 7 = cos 9 cos 6'' + sin 9 am9' cos(0 — 0')] and c is a posi- 
tive constant, and Gi^{r, r') is the 'radial' Green function 
satisfying 



2 ,/(L2_i/4) 2M\ 



(11) 



together with ingoing (outgoing) boundary conditions at 
the horizon (spatial infinity) and a causality condition. 
The radial Green function may be constructed from a 
pair of linearly-independent solutions (mJJJ, [defined in 
Eq. ^] and u"p [defined by zC - e+'"''- as 00]) 
of the homogeneous equation mn, via 



G,^(r,r') = -- 



W{l,uj) 



(12) 



where r*> = max(r*,r^), r*< = min(r*,r^) and the 
Wronskian W{l,oj) is 



Wil,L0)=K 
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(13) 



Note that the Wronskian is zero at QNM frequencies, 
i.e. W{l,ujin) = 0. Equivalently, the 'in' and 'up' QNM 
solutions are linearly-dependent. Hence the QNM fre- 
quencies correspond to poles of Gi^^ . 

At suitably 'late' times t — t' > r^+r'^, the contour of 



integration defining the inverse Laplace transform ( 10 ) 



may be deformed in the lower half-plane, as shown in 
Fig. g The QNM contribution G?,f ^ arises from apply- 
ing Cauchy's integral theorem to pick up the contribution 



from the residues of the poles at QN frequencies {w/„}. 
It may be written in the following way: 
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Here the 'reflection time T' is defined by 

T = t — t' — r^, — r[, 
the normalised wavefunctions M;n(r) are defined by 



iiinir) = uZ^Jr) X 



(14) 
(15) 

(16) 



and the so-called 'excitation factors' Bi, 

by 



Bln = 



A 



(+) 



dA 



2uj 



duj 



are defined 



(17) 



The 'QNM sum' expression, given by Eq. ( 14) and aris- 
ing from the residues of the poles (2) in FigTjS was con- 
sidered in some detail by Leaver [50]. It is not the only 
contribution to the Green function; observe from Fig. [3] 
that we must also take into account the parts of the con- 
tour (1), the 'high-frequency arcs', and (3), the 'branch- 
cut integral'. Parts (1) and (3) are associated with the 
'prompt response' (direct propagation) contribution, and 
'power-law tail' (back-scattering from potential) contri- 
butions to the Green function [40l [41] , whereas part (2) 
is associated with 'quasinormal ringing' related to the 
peak of the potential barrier, which is associated with 
the photon sphere. 



The QNM sum Eq. (14) seems ill-defined at 'early' 



times T < 0, due to the exponential divergence of the 
exp(— iwinT) factors in the sum over overtones n in 
Eq. ( [l4| . By contrast, at late times T > 0, it appears 
that the sum over n will converge, due to the exponential 
suppression of high overtones. Andersson [22 demon- 
strated the convergence of the sum over n for fixed r. 
Note however, that if the source does not have compact 
support this does not imply convergence of the integral 
in Eq. ([9]), since for any T there will always be points 
in the integration outside the radius of convergence of 
Eq. ( [l4| . Furthermore, Andersson showed that, subject 
to the above caveat, a combination of the QNMs and the 
contribution from the branch cut can form a complete de- 
scription of the wave evolution for T > (in other words, 
the prompt response [contour (1)] can be neglected in the 
regime T > 0). 

In principle, if we know the excitation factors Bim fre- 
quencies^ ioin and wavefunctions u;„(r) we can deduce 



from ( 14 ) the QNM response excited by a given perturba- 



tion. However, such calculations are relatively technical 
and have only been attempted in only a handful of stud- 
ies in the literature [20H24] . In this paper, we develop 
approximation methods in an attempt to investigate a 
physical question: what is the structure of the Green 
function close to the lightcone? 
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IV. LARGE-; ASYMPTOTICS OF QN MODES 

To analyse the large-Z asymptotics of the QNM sum 
Eq. (10), we need asymptotic approximations for (i) the 
QN frequencies, (ii) the QN wavef unctions and (iii) the 
QN excitation factors (QNEFs). These ingredients are 
obtained in this section. We set M = 1 in much of the 
following. 



Note that the ansatz has a geometric motivation; the 
phase factor in Eq. ( 18 1 closely resembles the (square 



root of the) right-hand side of the geodesic equation in 
the critical case b = be = V27 (where b — l/uj is the 
impact parameter for a geodesic approaching from spatial 
infinity), i.e. 



A. The expansion method 

In Paper I, we described the key steps in the expansion 
method. The first step is to introduce an ansatz for the 
wavefunction, 



uiuj{r) = exp 



lOJ I I 1 



1/2 



vi^{r). 



(18) 



T.4 



= 1 



1 - 



(19) 



Note also that the QN boundary conditions are automati- 
cally satisfied if vi^{r) is regular in both limits — > ±(X). 
Now, substitution of ansatz ( 18 ) into Eq. ([6| leads to the 
equation 



ifv' 



2icj 1 + 



1/2 



27^2 - ^ 27iw . 6 



-1/2 



2/3 



4r2 



u = 0, 



where here ' denotes differentiation with respect to r. The equation above looks more complicated than Eq. ([6|, but 
it has the advantage of being amenable to an expansion in inverse powers of L. We expand the frequency and the 
wavefunction in inverse powers of L as follows: 



J=-l 



(20) 



and 



vin{r) 



r 



n n oo 



i=l j=l 



exp(S'o„(r) +i iS'i„(r)-f...) 



(21) 



Here Wjn and aijn are complex coefficients, and Sj„{r) 
are smooth functions of r. Inserting the ansatz into 



Eq. (20) and collecting together like-powers of L leads 



to a system of equations. To solve the system, we im- 
pose a condition of regularity on the functions Sjn at 
r — 3M. This leads directly to the frequency expansion 
Eq. ([2]) and to the results quoted below. More detail on 
the expansion method is given in Paper I, and in Ap- 
pendix |A) 



exp iLUin 



B. QN wavefunctions 

To leading order in L, the QN wavefunctions are 
Ui„(r-)«C„exp(5o„(r))(l-3/r)" 

(l + 6/r)^/^(l-3/r)-l f-^dr 

(22) 

where C„ is a normalisation constant such that 
limr->.oo uiniT) = L The function SQn{r) is found to be 

Sonir) = lH2/x) + 2{n + l/2)ln ( ) (23) 

2 V .T + v3 / 



where 



a; = 1 



1/2 



(24) 
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We makes three notes in passing: (i) x varies between 
1 (at spatial infinity) and 2 (at the horizon), (ii) the 
constant of integration in finding Son{f) was chosen such 
that Son{r = 2) = 0; (iii) to leading order the radial 
wavefunctions are independent of the spin of the field. 
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The integral in Eq. ( 22 I can be found in closed form, 



exp 



(l + 6/r)i/2(i_3/r)-ll f'^dr 



(25) 



hence 



^*;n(r)«(l-3/r)"x-i/V'"^K-i)^ 



1 + a; 
4- 2.T 



4iuj 



i + vr 



2n+l 



(26) 



Inserting the lowest-order approximation for the fre- 
quency. Win = (27)~^/^(L — iN), the QN wavefunction 
may be rewritten in the following form: 

uin{r) « U{rMr)r exp (^L7^(r)/^/27) (27) 
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FIG. 4. Excitation factors Bin of fundamental (n — 0) modes 
of scalar field. The plot sliows the real and imaginary parts 
of the excitation factors Bto for modes Z = 2 ... 10. The blue 
dotted line and markers show the asymptotic approximation 
[Eq. m}]. The red crosses show the numerically-determined 
values (see e.g. [24) V 



where 



U{r) = ^-i/2e-K--i)V(4y27) 



1 + x 
4- 2a; 



2/\/27 



1 + Va' 



7^(0=41n(i^)-l.(--l)^ 



p{r) 




(28) 
(29) 



(30) 



where the 'geometric constants' are 



C = 3-V27 + 41n2-61n(2 + V3) 

« -7.325311084, 

./27e';/V27 

B = = « 0.19182703317, 

(2 + V3)\/^ 

K = ^— « 0.92482482643. 

(2 + 73)2 



(32) 
(33) 

(34) 



Excitation factors 



Fig. |4] shows the 'fundamental mode' [n = 0) excita- 
tion factors for the scalar field, in the complex plane. 
It compares the lowest-order large-/ approximation given 
in Eq. (31 1 [blue] with numerically-determined excitation 
factors [red] , for modes / = 2 . . . 10. The agreement is 
good, and improves as I is increased, as expected. 



The excitation factors defined in Eq. (17) are some- 
what more difficult to obtain, since we must compute 
the derivative of the ingoing coefficient with respect 
to w, and evaluate at the QNM frequency. To leading or- 
der in L, it is straightforward (but somewhat tedious) to 
compute the excitation factors by combining the expan- 
sion method with standard WKB techniques. We present 
the full analysis in Appendix |Aj and here merely quote 
the key result: 



•l/2^-l/2_gg2iCL/y27 



nl 



(31) 



V. STRUCTURE OF THE GREEN FUNCTION 

In this section we employ the large-Z approximations of 
the previous section to investigate the singular structure 
of the Green function; in other words, the form of the 
GF close to the lightcone. Let us begin by rewriting the 
QNM sum (fT4| as 



^rot (X,X ) — 



2 

rr' 



Re'^gi{x,x'), 

1=0 



(35) 
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where 



gi{x,x') = ^{21 + l)Pi{coaj)BinUin{r)uin{r')e'' 



Tl = 



iuiinT 



(36) 



Note that the series (35) is divergent, in the sense that 



limj_5.oo l^^il 7^ 0. In the large-/ hmit (and for fixed 
points x') the phase of Qi varies linearly with I, i.e. 
lim/_j.oo ai'g(C//+i/CJi) = const. In Sec. VB we show that, 
in general, well-defined values can be extracted from the 
series ( 35 1 using an appropriate summation method j42j . 



However, the summation method does not yield a finite 
result if the terms Qi combine coherently as / — >■ oo. In 
this case, the series representation of the Green function 
is singular. 



A. Singular structure and the light cone 

As argued above, the QN mode sum for the Green func- 
tion is singular if the coherent phase condition is satisfied: 



lim a.T:g{Qi+i/Qi) = 27rfc, 



fc e Z. 



(37) 



From physical considerations, we expect the Green func- 
tion to be singular on the light cone (i.e. if the spacetime 
points X and x' are connected by a null geodesic). Hence 
we expect to find that the phase of Qi is related to the 
null geodesies on the black hole spacetime. Below, we 
use the asymptotic results of Sec.|TV]to demonstrate the 
link explicitly for the Schwarzschild spacetime. 



Inserting asymptotic results (27 1 and (31) into (35) wc 



obtain that, at leading order \\\ L = 1 + Tj2 
lim Gi = 2BU{r)U{r')e-^/^^^\iLy/^ 

X ^^L[(2C+n(r)+^Z^r')-T)/V27]p^_^J^^^^■^ 

°o [(-iKp(r)p(r')ie-^/^)" + C'(L"-i 

n=0 



(38) 



If we are prepared to neglect all subdominant terms in L 
then the sum over overtones n may be performed. 



n=0 



exp {—iKp[r)p[r')Le 



-T/\/Tf 



(39) 



Let us now assume that 7 7^ 0, tt, and employ the ap- 
proximation 



lim Pi (cos 7) 

/7»1 



/ 1 \ 1/2 

( ) f e-"/4g»i7 + ^^^/A^-^L^'\ 

\ 2itL sin 7 / V / 



(40) 



With results (381, (39) and (40), the coherent phase con 



dition ( 37 1 is equivalent to the implicit equation 



T 



27 {2nk ± 7) + 2C + 7^(r) + 7^(r') 
- V27Kp(r)p(r')e-^/^, 



(41) 



where k is an integer. Quantities T, K,TZ(r), p(r) were 
defined in Eqs. p^, (IMl), (ISOl), respectively 



The quantity 2TTk ± 7 (where < 7 < is the angle 
between the spacetime points) may be interpreted as the 
angle traversed by a null geodesic in an orbit around the 
black hole (see Fig. [T]). 

We expect the solutions of Eq. (41), i.e. T — 



Tl^'^^^\r,r' to be approximations to the times at 
which the Green function is singular; and on physical 
grounds we expect the Green function to be singular if 
(and only if) x and x' are connected by a null geodesic. 
Hence we may test our argument by comparing j'^Q^^'* 

against the 'null geodesic times', found by inte- 

grating the Schwarzschild orbital equations. 

Table |l] compares the 'QNM times' with 'geodesic 
times', for the special case of 7 = and r = r' . In 
the first column is the time T^^^eo) — /^^ _ where 
At is the coordinate time it takes for a null geodesic to 
start at a point, orbit once around the black hole, and re- 
turn to the same point (i.e. r = r', 7 = 0). In the second 



column is the time y^Q^^^) found by solving Eq. (41 1, us- 



ing a simple iterative method. We find excellent, though 
not precise, agreement. The 'QNM time' approaches the 
'geodesic time' as the number of orbits goes to infinity 
(i.e. k — >■ 00). 





y(geo) 


y(QNM) 


y(goo) _ y(QNM) 


r = 4 


25.6449315 


25.6449191 


0.00001246 


r = 10 


20.7019188 


20.7017749 


0.00014394 


r = 100 


18.1280305 


18.1285177 


-0.00048718 



TABLE I. Singularity Times. Here T = At- 2r, where At is 
the time taken for a null geodesic starting at radius r to orbit 
the black hole once and return to the starti ng r adius. y(QNM) 
is found by solving the implicit equation (41 1, and T^sco) jg 



found by solving the geodesic equations numerically. 

Fig. [5] shows a cross section [in the plane = 7r/2] of 
the wavefront at t — ti associated with an initial event 
at x' = 8M, y' = t' = (see also Fig. [5]). Here we have 
chosen ti « 41. 8M, which is the time is takes for the 
wavefront to orbit once and return to a; = 8M. This is 
sufficient time for two caustics to form. The dotted blue 
line shows the prediction of the position of the wavefront 



found from inserting the QNM asymptotics of Sec. IV 
into the coherent phase condition Eq. (37) [more explic- 
itly, for a given fc, ±,7 we solved (41) numerically and 
rearranged (15) to find r.^, and hence r]. The plot shows 



that the wavefront prediction is most accurate for large 
T, in other words, close to the unstable orbit at r = 3M. 
It is natural to presume that the (exact) QN mode sum 



( 10 ) is in fact singular precisely on the light cone (and not 
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Equatorial cross section of light cone from point (x=8, y=0) at time t = 41 .837 M 



30 
20 




-30 



-10 



10 
x/M 



20 30 



50 



FIG. 5. Wavefront cross section. This plot shows a 'snap- 
shot' in the equatorial plane of a wavefront [red solid line] 
emanating from x' — 8M, y' — 0, t' — 0, taken at coordinate 
time t — ti ^ 41. 8M. It corresponds to a cross section of the 
light cone shown in Fig.[l] Note that the wavefront intersects 
itself at two antipodal points (caustics). The blue dotted line 
shows the position of the wavefront predicted using the large-? 
asymptotics of the QNMs [Eq. ^]. See also Fig.[7| 



just close to it). That is, we presume that the difference 
between the QNM time determined from Eq. ( [41] ) and 
the geodesic time must arise because we have neglected 
important sub-dominant terms in the expansions in L. 
In particular, performing the sum over n in Eq. (39) is a 



non-rigorous step that should be viewed with some sus- 
picion. It may be possible to improve Eq. (41 ) somewhat 
by extending the approximations (27) and (31) to higher 



order in L. 



B. Caustics and the four-fold singular structure 

In the previous subsection we demonstrated the link 
between the light cone and the large-Z asymptotics of 
the terms in the QN mode sum ( 14 ) . Here we make 



further use of the large-/ asymptotics to investigate the 
form of the scalar-field Green function close to the light 
cone. We show that the singular part of the Green func- 
tion undergoes a transition each time a caustic is en- 
countered (where we use the term caustic to refer to a 
self- intersection of the light cone) . 

Let us start by inserting the large-/ asymptotics of 
Sec. |IV| i nto the QNM decomposition of the Green func- 
tion ( 14 ) . We obtain the expression 



OO 

Re^(zi)i/'P,(cos7)e-*^*('-'''''T'^) (42) 



1=0 



where 



X{r, r', T) = ^ ^ (43) 

rr' 



and 



$(r,r',7,T) = 



^-7^(r)-7^(r')-2C 



27 



Kp{r)p{r')e 



(44) 



Here Grot^ ^ quantity which should have the same 
singular structure as Cji^^'^ and (we assert) Gret- 

We now follow the steps taken in Sec. VD of ^331 , where 
a similar analysis was performed for the Nariai spacetime 
[33] ■ Related methods are used in, e.g., seismology [33] ■ 
First, we make use of the Poisson sum formula to convert 



the infinite sum (42) into an integral, viz., 



G, 



(sing) 
ret 



X 



+ 00 



(-l)'^Re / dLe^'''''^(iLfl'^e-''^'^Pi^_^i-2{co&l)=Xy^Z, 



L=0 



m=0 



(45) 



where 



and 



I,n = Re / dL{iLfl^e-'^'^Rn 
Jo 



(46) 



m even. 



(-l)('"+i)/2 (Q^+\/2(cos7)e^(™+i)'^^ + Q^l\/2(cos7)e-^(™+^)''-^) , m odd 



(47) 
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with 



QW(cos7) = I (pm(cos7) ± ^g^(cos7)^ 



(48) 



where Q^{-) is the Legendre function of the second kind. the following simple 'exponential' approximation: 



The integrals ( 46 ) can be approximated if we make use 

functions 



of asymptotic approximations for the 
which are valid in the large-Z limit. We will make use of 



)(±) 

-L-l/2 



(cos 7) 



1 



27rL sin 7 



1/2 



„±i7r/4 =FiL7 



(49) 



This approximation is valid in the regime L7 ^ 1 (and 
L\'K — 7I ^ 1); in other words, it is inappropriate for use 
if the spacetime points are separated by a small angle 
7 ~ 0, for example, close to a caustic. 



Applying approximation (49) we obtain 



1 



27r sin 7 



1/2 



Re 



dL 



L=0 



(-l)™/2 



_L(m7r+7 — 



-iL{ 



|^jgiL([m+l]7r-7-*) _|_ g-iL([m+l]7T-7+*) j 



m even, 
m odd, 



(50) 



for L ^ 1. The integral will be singular if the phase 
factor in either term in ( 50 ) is zero; unsurprisingly, this 



is equivalent to the coherent phase condition we inves- 
tigated in Sec. VA The QN mode sum (10) is only ex- 



pected to be convergent for T > (hence $ > 0) and so 
for simplicity let us neglect the latter terms in ( 50 ) . Now 



we may regularize the integral by introducing a small 
positive imaginary part in the exponent [introducing a 
small e into the exponent of ( 50 ) is equivalent to impos- 
ing a smooth high-Z cutoff upon the integral], and use the 
identity 



lim 

f-fO+ 



-+TTSiq), 

q 



(51) 



to obtain 



2 sin 7 



1/2 C (— 1)™/^ (5(m7r + 7 — $), m even. 



(_l)(™-Hi)/2 



7r[$ - (m + l)7r + 7] ' 



m odd. 



(52) 

From the above expression, we see that the sum is singu- 
lar when $ = 27rfc ± 7 (where k is an integer) , as antic- 
ipated in Sec. VA [see Eq. (41)]. However, the form of 



the contribution I,„ to the Green function depends cru- 
cially on whether m is even or odd. The 'shape' of the 
singularity alternates between a delta distribution and a 
singularity with antisymmetric 'wings'. The mth wave 
may be associated with the mth orbiting null geodesic. 
Now, m has a clear geometrical interpretation: it is the 
number of caustics (antipodal points) through which the 
corresponding geodesic has passed. Equation ( 52 ) im- 



plies that the Green function changes form each time the 
wave front intersects itself (i.e. passes through a caustic) . 

The simple argument above gives us some insight into 
the 'four-fold' singular structure of the Green function 
|32j . However, the approximation (49) leading to (52) 



breaks down when 7 is close to or tt. For example, the 
approximation is not valid close to a caustic. To improve 
our analysis, we may employ the uniform asymptotics 
estabhshed by Olver [45] . 



2i?i/2(cos7) 



7 
sin 7 



1/2 



H^ijL), (53) 



where Hq{-) — Jo(-) ± *^o(') are Hankel functions of the 
first (-I-) and second (— ) kinds. The asymptotic approxi- 
mation (|53| is valid for / 1 for all angles 7 < vr. Using 
([53| in (|46| leads to 



1/2 



(_1)W2 [i/+(7L)e*-^("'^-^*) +i/o-(7i)e-^^(™'^+*)] 



TO even, 



Re / dL{iLy/^ 

4sm7y Jo I (_l)(m+l)/2 [^-(^_^)g^L([m+l].-<I.) ^^+(^^)g-^L([m+l].+$)] ^ ^Ojj^ 



(54) 



Such integrals can be found with the aid of results in standard tables books, for example, Gradshteyn & 
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Ryzhik [15]. Once again, for simplicity let us neglect 
the latter terms (which are singular for T < 0). Con- 
sider now the even-m case; the first term is singular when 
$ = rmr + 7. If <I> < rmr + 7, we may integrate along 
the positive imaginary axis by making the substitution 
L = iv, and using the identity HQ{i^v) — 2Ko['^v) / [iir). 
The integral is purely-imaginary, and hence upon taking 



the real part we get zero. Conversely, if $ > rmr -f 7, 
then we may integrate along the negative imaginary axis 
by making the substitution L — —iv and using the iden- 
tities H+{-i-iv) = 2I0M - H^f^-i-yv) = 2/0 (7i^) + 
2Ki^{-^v)/{i'K). Taking the real part picks out only the 
latter term, and the integral may be found with the aid 
of Eq. 6.621(3) in [151. Proceeding in this way, we obtain 



1 



7 



2 \sin7 



1/2 



/ \ 1/2 

— (5($- [m7r + 7]), $ < m-Tr + 7, 
(-1)W2 J V / 

2^^i(i,i;2;|E^), $>m. + 7, 



2[<i> — TOvr -I- 7 



3/2- 



m even, 



-2E (27/[(m + l)7r + 7 - $]) 



(_l)(™+i)/2 J + l)7r - 7 - $][(m + l)7r + 7 - $]i/2 ' 



$ < (to -I- 1)tt — 7, 



2[$ - (m + l)7r + 7]3/2 



F.{lh^; %XXZ7, )^ $>(TO-Hlh-7, 



(55) 



TO odd. 



Here 2F1 is the Gauss hypergeometric function, and E 
is the complete elliptic integral of the second kind. Note 
that the delta function deduced in Eq. (52) has been 
included in the top line of Eq. ( 55 1 . Eq. ( 55 ) is valid 



in the range < 7 < 7r/2; but extending the result to 
7r/2 < 7 < TT is straightforward. 



Numerical visualisations 



Figure [6] shows the asymptotic approximation (55), for 
fixed points separated by angle 7 = 7r/24 at radius r = 
r' = 8Af as a function of reflection time T = t—t'—r^—rl. 
The contributions from Ii and Z2 are shown separately. 
The GF is singular at T = Ti, T2, . . ., where Ti and T2 
correspond to the time taken for a null geodesic to pass 
through an angle 27r — 7 and 2tt + 7, respectively. The 
former geodesic (1) passes through one caustic, and the 
latter (2) through two. Near Ti the GF is approximately 
antisymmetric. At T2 the GF has delta-function support. 
Beyond T > T2 the contributions from Xi and X2 cancel 
precisely. 

The structure of the GF in the xy plane (i.e. 9 = 7r/2) 
is illustrated in Fig.[7j It shows a snapshot att = ti of the 
wavefront [solid lines] and GF [shading] associated with 
an event at x' = 8M, y' = t' = (see also Fig. [5]). The 
first and second caustics are visible (at x « — 24.36M 
and at a; = 8M). Observe that different parts of the 
wavefront have passed through zero (pink line, outer), 
one (green) and two (purple) caustics. The plot illus- 
trates how the GF changes depending on the number 
of caustics crossed. On the outer wavefront (0 caustics, 
pink), the GF has delta- function support. Around the 
intermediate wavefront (1 caustic, green), the GF is an- 
tisymmetric (positive on the exterior, negative on the 
interior). On the inner wavefront (2 caustics, purple). 



r = r' = 8M 
Y = It/24 


j 
i 

\ 




\ 


\ 


(2) 



T = t - 1' - r. - r.' 



FIG. 6. Asymptotic approximation to the singular structure of 
the Green function. This plot illustrates the GF as a function 
oi T — t ~ t' — r, — ri , for fixed spatial points at r = r' = 8M 
separated by a small angle, 7 — 7r/24. The black dashed 
line was computed using approximation to the Green func- 
tion given by Eq. (45 1 and ( |55[ ). The red curve (1) shows the 
contribution from Ii, and the blue curve (2) shows the con- 
tribution from I2. The 'spike' at T ~ 22.13 indicates a (-ve) 
delta function. In the regime T > 22.13, the terms Xi and X2 
cancel precisely. 



the GF again has delta-function support, but with the 
opposite (-) sign. 



VI. DISCUSSION AND CONCLUSION 

In this paper we have used the expansion method 
of [13] to investigate the large-/ asymptotics of the 
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ansatz f3IP for the 'direct' part of the GF, i.e. 



FIG. 7. Illustration of the singular structure of the Green 
function. This plot shows a 'snapshot' of the scalar-field 
Green function (GF) in the xy plane a,t t — ti associated 
with the initial point at x' = 8M, y' = 0,t' = 0. Note that 
ti ~ 41.837Af is the coordinate time it takes for part of the 
initial wavefront to orbit the BH once and return to x = x' , 
y = 0. On the outermost wavefront [pink line], which lies 
outside the first caustic point at a; « — 24.36M, y = 0, the 
GF is proportional to a positive delta function [Eq. |50|]. On 
the innermost wavefront [blue line], inside the second caus- 
tic point (at X — 8M,y = 0), the GF is proportional to a 
negative delta function. Across the intermediate wavefront 
[green line] which joins the two caustics, the GF has antisym- 
metric 'wings' approximated by 1/a (where a is the Synge 
worldfunction) . The dotted black line shows the position of 
the wavefront estimated from the lowest-order asymptotics of 
the QNM sum (see also Fig. [5|. The red (positive) and blue 
(negative) shadings indicate the value o f th e GF (computed 
via the asymptotic approximations [Eq. (55l]). The coloured 



disk indicates the area of spacetime in which the sum over 
QNM overtones [Eq. (|36|l] is convergent. 



Schwarzschild QN modes and their excitation factors. 
We used the asymptotic results to investigate the QNM 
contribution to the scalar-field retarded Green function. 
This led to insight into the singular structure of the GF 
near the lightcone. We showed that the form of the GF 
changes every time a caustic is encountered, and that 
it undergoes a four-fold repeating pattern [see Eq. ( 50 1 , 
([52]), ([55]) and Fig.[7]. 

A four-fold singular structure for a 4D spherically- 
symmetric spacetime (like the Schwarzschild black hole) 
has been anticipated (e.g. [351 however, we believe 
this is the first time it has been demonstrated for a black 
hole spacetime using a QN mode representation. As was 
discussed in Sec. VE of Ref. ^3], an alternative way to 
make sense of the four-fold pattern is with a Hadamard 



'-'ret {^y^ ) 



lim Re 



Ai/2 



Re 



I 



(56) 



Here a — a{x,x') is the Synge world function [48 (half 
the squared interval along the geodesic joining x to x'), 
and A — A{x,x') is the van Vleck determinant, a mea- 
sure of the degree of focussing of neighbouring geodesies 
[5 5) . The determinant is singular at the caustic. In pass- 
ing through a caustic on the Schwarzschild spacetime, A 
changes sign. If we presume that this leads to an addi- 
tional factor of —i in the square bracket of ( 56 1 each time 



a caustic is traversed, we recover once more the four-fold 
pattern S{a), l/(7rcr), —S{cr), — l/(7rcr), S{(t), etc. 

The effect of caustics upon wave propagation has been 
studied with a variety of methods (e.g. WKB methods, 
path integrals, functional integration, etc.) in a range 
of disciplines including optics [5D], acoustics [5T], seis- 
mology [44] , symplectic geometry [52] , and quantum me- 
chanics [531 \SM ■ The formation of caustics in black hole 
spacetimes has also been studied extensively, in the con- 
text of strong gravitational lensing [5SH5S] . Some further 
investigation of the effect of caustics on the details of 
wave propagation (i.e. on the Green function) in curved 
spacetimes now seems to be warranted. In particular, 
analysis of the singular structure of the Green function 
near a rotating (Kerr) black hole, an axisymmetric space- 
time, would be of great interest For this challenge, 
the utility of the QN expansion method appears rather 
limited [T5]. Instead, let us focus (for now) on what more 
can be done in Schwarzschild with our methods. 

In Paper I [T3] we introduced the QNM expansion 
method and applied it to investigate the QN frequency 
spectrum. In this work - Paper II, let's say - we have 
explored the geometric interpretation, using first-order 
(in L) expansions to bring qualitative understanding of 
the singular structure of the GF. It seems there is scope 
for a future work (i.e. Paper III) in which higher-order 
(analytic) expansions are combined with numerical meth- 
ods [ID. 24 in the small-/ regime, to achieve quantitative 
calculations of QNM sums [such as Eq. (14)]. 

Let us suggest two areas where this work may find 
further applications. The first is in the calculation of the 
gravitational self-force [331|S3], which is required to (e.g.) 
model extreme mass-ratio black hole inspirals [57], test 
the cosmic censorship hypothesis [58] , and calibrate effec- 
tive one-body theories [55]. With complete (i.e. global) 
knowledge of the GF, it is straightforward (in principle) 
to compute the self-force; yet the former is not trivial to 
obtain. In the 'method of matched expansions' [60j one 
seeks to match a 'quasilocal' expansion of the GF [61| . 
valid inside a normal neighbourhood (approximately, up 
to the first caustic), to a 'distant-past' expression. The 
QNM expansion makes up part of the latter; however, 
it must also be augmented by an accurate calculation of 
the branch cut integral. This remains to be achieved. 
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A second possible application is to the theory of Com- 
plex Angular Momentum (CAM) on black hole space- 
times [5^H55] . It is well-established (see e.g. [BB]) that 
the black hole photon orbit creates regular interference 
features in scattering cross sections, that is, 'spiral scat- 
tering' oscillations at intermediate angles and a 'glory' in 
the backward direction. In [S7] it was shown that such 
oscillations are intimately linked to the Regge poles of 
CAM theory. Regge poles are the complex-? versions of 
QNMs. Very recently, Decanini and Folacci [68] have ex- 
tended the expansion method to obtain high-frequency 
asymptotics for the Regge pole residues (which, in CAM 
theory, play an equivalent role to the QNM excitation 
factors Bin)- We anticipate that their work will lead to 
improved analytic estimates for spiral/glory scattering. 

In conclusion, we hope that our study of the large-Z 
asymptotics of QN modes has shed some light on the 
global structure of the Green function on a black hole 
spacetime, and we anticipate further illumination from 
work to come. 
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Appendix A: QNM excitation factors Bin via WKB 
analysis 

In this section we derive a leading-order (in L) approx- 
imation to the excitation factors, defined in Eq. (17 1. 



The key ingredient in the calculation is the derivative 
dA[J /du!, taken at the QNM frequency. To find this, we 
must perturb the frequency slightly away from its QNM 
value, uJin = (jJiri + e, (where e is small) to determine 
the first-order change in a\J . Note that we use overbar 
notation uj to denote the precise QNM value. In other 
words we perturb only tzjq, the 0{L^) term in expansion 
pO|, so that 



Too = -iiV/V27 + e 



(Al) 



This results in a breakdown of continuity at r = Tc = 3 in 
expansion (21 1. However, we may find regular solutions 
in a 'interior' region close to r = 3, and match these 
onto 'exterior' solutions (in regimes r > 3 -I- e and r < 

3 — e), to obtain A\J = 0{e) to first order. In short, 
we combine the new expansion method with a standard 



WKB approach to derive the leading-order result for the 
excitation factors. 



1. Interior solution 

To find an 'interior solution' valid near r = 3, let us 
begin with the radial equation ([6]) and make the ansatz 
uiuj{r) — so the radial equation becomes 

^ • U{r)^P = 0, 



where 



C/(r) = /- 



L2 - 1/4 



M 



(A2) 



(A3) 



Now substitute in the expansion (20) for a;, and make 

(A4) 



the change of variables from r to z via 



so that 



U{r) = {2V3 mo + z^/3)L + (l^/^^ . (A5) 



Then, at leading order (L^) the differential equation (A3) 
implies 



(2V27ujo + z^^tP = 0. 



(A6) 



Now substitute in the perturbed frequency ( Al I to obtain 

^ + + 2y27e) V = 0, (A7) 

which has independent solutions 

V'l = Ai+r,((-l+*)2), i^2^ Dn+„{{l-t)z) (AS) 

where r] = ^/27ie and Da{ ) is a parabolic cylinder func- 
tion [47]. We will use "ipi, since it turns out to have the 
correct 'ingoing' behaviour towards z — >■ —00. The large- 
z asymptotics of the solution ipi are straightforward to 
obtain. As z — > —00 we have 

^ 2("+'')/2e^*^("+'')/4|z|"+''e+'^'/^ (A9) 

and as z — > +00 we have 

_ (27r)i/^ 

X 2-(i+'?+")/2|^|-('7+n+i)e— V2. (Aio) 
Hence to lowest order in e these are simply 

- 2"/^e-^"'/^|z^e+^^"/^ z^-oo, (All) 

and 

"/2 



g't7r(n+?7)g— 3i7r(l+ri+jj)/4 



e V27 z r(n + l)(2^)i/2e-3-(«+i)/4 



(A12) 
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2. Exterior solutions 



where 



In this section we construct solutions vahd in the ex- 
terior regions, i.e. on either side of the interior around 
r — 3. Let us define solutions 

u=^(r) =exp ^±«^^ {l + 6/r)^^'^{l-3/r)dr^ w=^(r) 

' (A13) 
where satisfy Eq. and are solutions of 



dr 



r 

,2 r2 



2 / 6 

± 2ia; I 1 



1/2 



27w^ - 27ioj f 6 
^ 1 + - 



-1/2 



dv"^ 
dr 

1 2/3 

^3 



Z'(r) =\/27/ (r-3)V7-(?- + 6) (A22) 

Z{r) =ln|r-3| 

-In (3 + 2r + ^/3r{r + 6) ) +ln^ 

(A23) 



Let us now make the change of variables ( A4 1 and exam- 
ine close to z = 0: 



S+ w In 



=0 



(4/3)i/4(C/12)"+i/2 
Zwln \z\/{2Vf7L) 



(A14) 



Note that, around the critical orbit r = 3, we have 



(A24) 
(A25) 



The asymptotics of the and u functions defined in 
dAlSl) are 



exp(±iio {l + 6/ry/^{l-3/r)dr. 



exp(±izV2 + C'(i"^)). 



(A15) 



As in Eq. (21), we make the ansatz 



= [(1 - 3/r)" + O(L-i)] exp(5o± (r)+L-i5± + . . . ). 

(A16) 



After plugging in to ( A14| and matching like powers of 
L we obtain at order the following: 



,4^27/ \2V27L 



27/ V2V27L 



(A27) 



±2i\l + - 
r 



1/2 



6m 



1- j5o"n±^|l 



1/2 



^ 27i / ^ , 54t^i 



0. (A17) 



At the QNM frequency = wq = —iN/\/ 27 the solu- 
tion is SQ^{r) — SQ^{r), where 



1 



H^r/2) - ln(3 + 2r + v/3r(r + 6) ) 



(A18) 

Here for convenience we've defined a (seemingly ubiqui- 
tous) constant 



^i^ = (2 + ^/3f = 7 + 4^/3. 
2 — v3 



(A19) 



Now let us perturb about the QNM frequency, employ- 



ing ( Al ), to obtain 



SL (r) = S+„ (r) + V27ieZ{r) ( A20) 

Sonir) = SUr) - [(2n + 1) + ^/2riie] Z{t) (A21) 



3. Matching procedure 



Let us now define three solutions, 



W0 = £'„+V27ie((-l+«)^) 



Uy = Bin" + BoutU^ 



(A28) 
(A29) 
(A30) 



Now we match uq to u< in the region z < 0, and match uq 
to u> in the region z > 0. Applying asymptotic results 
(lAllL (|A12|), (|A26| and (|A27|) we obtain 



Cin ^22"2("+i)/27^"/2 r^/^j g-iW4 (A31) 



=(-i)"an 

Bin=er(n-f 1)2-"L-" 



-N 



1/2 



X e 



-iTrn/2 



2V27 



2N 



(A32) 
(A33) 



(A34) 
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4. Ingoing and outgoing coefficients 



Excitation factors 



In order to find the ingoing and outgoing coefficients 
J and A''i^\ we first need the foilowing 'phase factors': 

0^/;=t(l + 6/r)i/2(i _ 3/r)/-idr) 

(A35) 



exp 



ai = 



Pi 



7i 



exp(— iwr*) 
exp (^iuj[& - \/27 + 81n2 - 3hr^] 

exp (»c./;r3°°(l + 6/r)V^(l - 3/r)/-Mr) 
exp(+ia;r*) 

exp(^iw[3-\/27 + 41n2-31n,f]) (A36) 
exp + 6/r)i/2(i _ 3/r)/-Mr) 



exp(— zcjr.^) 



= l//3i 
a2 ~ iim e 



TV 



-N 



h= lim 6^0+^21/2(^1/72) 

r— ^00 \ / 

72 = lim e^'^" = 2^/2 (2,^1/2 

r— f 00 V 

Then, to leading order in L, 

^(+) ^ /3l^2 -Bout 

(-) 7172 -Sin 



(A37) 
(A38) 

(A39) 

(A40) 



N 



(A41) 



A 



aiaa Cin 



X e 



-iw[9+2y27+121n2-61n5]2l/2 |^2^1/2^ 



2V27 

N 



2N 



(A42) 



The excitation factors Bin are defined in Eq. (17 1 



Combining results (A42) and (A41) we obtain the 



leading-order approximation 



ri-1/2 



exp(2iw;„C) 



n! iu;iJ[V27L]) 



Stt 



216 



AT 



(A43) 



with error 0{U^ '^/2) and where the constant Q is defined 
in (32 1. This can be rewritten as Eq. (31 1, the result given 



in the main text. 



Expression (31 ) fits the large-Z numerical data f55 well, 
for n = 0,1,2. It is considerably more effort to obtain 
the higher order corrections and we shall present these 
together with efficient numerical techniques for their cal- 
culation in a subsequent publication. 
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